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A non-relativistic scalar field coupled minimally to electromagnetism supports 
in the presence of a homogeneous background electric charge density the existence 



of smooth, finite-energy topologically stable flux vortices. The static properties 



^ ■ of such vortices are studied numerically in the context of a two parameter model 
describing this system as a special case. It is shown that the electrostatic and the 
mexican hat potential terms of the energy are each enough to ensure the existence 
of vortex solutions. The interaction potential of two minimal vortices is obtained 
for various values of the parameters. It is proven analytically that a free isolated 
vortex with topological charge N ^ is spontaneously pinned, while in the presence 
of an external force it moves at a calculable speed and in a direction (N/\N\) 90° 
relative to it. In a homogeneous external current J the vortex velocity is V = —J. 
Other theories with the same vortex behaviour are briefly discussed. 
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1. Introduction 



We study in detail some aspects of the physics of flux vortices in systems de- 
scribed by a non-relativistic complex scalar field minimally coupled to electromag- 
netism and in the presence of a background homogeneous electric charge density. 
The model was introduced and the essential features of its vortex dynamics were 
discussed in a recent publication [1] . Here we elaborate upon and extend the results 
of reference [1]. 

Our motivation to analyze the behaviour of topological solitons in the above 
model is twofold. First, it is actually a whole class of models potentially relevant 
to the theoretical treatment of several physical systems. In the absence of any 
scalar potential one deals with a plasma with infinitely-massive positive charges 
frozen in a configuration of constant charge density. With a quartic potential for 
the scalar field the theory can be viewed either as the natural coupling to the elec- 

[2] 

tromagnetic field of the dynamical Gross-Pitaevskii equation of superfluidity or 
as one possible dynamical extension of the static Ginzburg-Landau model of su- 
perconductivity, reminiscent in particular of the large friction limit of the system 
considered in reference [3]. Second, in the context of the above models one can 
obtain definite theoretical predictions about the gross features of the motion of 
vortices in two spatial dimensions, without any approximation and insensitive to 
the fine details of the Hamiltonian. This might prove a reasonable starting point 
for the understanding of the dynamics of an isolated Abrikosov vortex in thin su- 
perconducting films, a subject of considerable interest in the context of ordinary 
as well as high-T c superconductors 131 W W 161 . 

In section II we introduce the model 111 . Since the bulk of our results refer to 
two dimensional objects we directly present the model in two spatial dimensions. 
Section III contains our numerical study of axially symmetric vortex solutions for 
various values of the topological charge. We explicitly demonstrate the almost 
identical roles of the mexican hat potential term and the electrostatic term in the 
energy density of the model. We show that either one of these two terms is enough 
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to guarantee the existence of non-trivial vortex solutions. The parameter space is 
divided into two regions according to whether the ratio of the energy of the doubly- 
charged vortex over the energy of two well-separated single vortices is smaller or 
larger than one. Finally the two-vortex system is discussed and their interaction 
potential is obtained numerically. Contrary to the relativistic model the inter- 
action energy is not always a monotonous function of the distance between the two 
vortices. In a wide region of the parameter space there is a bump in the interaction 
potential of two vortices, which leads to the possibility of having doubly-charged 
dumb-bell shaped static solutions. The dynamics of vortices is discussed in section 
IV. The canonical structure of the model is derived and it is shown that the unam- 
biguous linear momentum is essentially the first moment of the topological density. 
The momentum conservation law in the absence of any external fields is shown to 
imply the spontaneous pinning of isolated vortices, while the application of an ex- 
ternal force results in the motion of the vortex at a constant calculable speed in a 
direction (N/\N\) 90° relative to it' 11 , features encountered already in the motion 
of magnetic bubbles in the ferromagnetic continuum' 91 1101 . Finally, in the discusion 
section V we show that necessary condition for the manifestation of the above 
Hall-behaviour in the dynamics of topological solitons in a given field theory is 
that the momentum part of its symmetry algebra admits central extension. Other 
interesting models sharing this property are commented upon. 



2. The model - General features 



We will be dealing with a non-relativistic complex scalar field \1/ (the conden- 
sate) minimally coupled with coupling q to the electromagnetic potential (Aq, A{). 
To make the model physically sensible and mathematically consistent it is nec- 
essary to introduce a background (positive-ion) charge density to neutralize the 
system. For simplicity we take it to be constant and homogeneous throughout. 
We concentrate on the physics of infinitely long straight vortices i.e. on field con- 
figurations uniform along the third spatial direction and define the model directly 
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in two space dimensions by the lagrangian 



2 

£ = ^[tt*P t * - c.c] + 7 ^oA)- ^-|^^| 2 + ^(E 2 -fi 2 ) - V(\*\) (2.1) 
z 2m o7T 



with 5 = etjdiAj , ^ = - diA , = (d t + iqA )V , = (9< - 

i|Ai)\l/ . 7, m and g are parameters, c is the speed of light and the spatial indices 
i, j range from 1 to 2. As it will become clear, although the specific form of the 
potential V changes the details of the profile of the vortex solutions, it does not 
affect their dynamical behaviour. For the discussion of vortex dynamics V could 
even be absent but in order to make contact with a variety of models of physical 
interest we will allow for a mexican hat phenomenological potential 



V(\V\) = I g (W - ^l) 2 (2.2) 



with quartic self-coupling g . Rescale fields and coordinates according to 



7« . . W_ . 

X 7 ^ Ob ii L ^ r» ~C 



*^* * A^^A,, a,^ C -^A, (2.3) 



10. 



to obtain' 11 , in terms of the dimensionless quantities t, Xi, Aq and A{ used from 
now-on, the lagrangian 

L = h^idt-Ao)* + c.c.) + A - I \D,$\ 2 + \{\v 2 - B 2 ) - \j 

2 2 2 A o 

(2.4) 

with B = eijdiAj , E{ = —dfAi — OiAq and D{ = di — iAi . We keep the same 
symbols to simplify our notation. 
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Classically the model depends on the two free dimensionless parameters k and 
A defined by 

9 9 3 4 

2 gm c m c 

The conversion to physical units for the lengths and the characteristic times is given 
for the specific system at hand by the rescaling formulas (2.3). An overall factor 
g^O^/n 2 was dropped from £ since it does not enter the equations of motion. It 
plays though the role of 1/h in the quantum theory and determines the necessary 
condition for the validity of our semiclassical approximation 



K 2 



9^ 



4 - (2-6) 





In this limit we expect the quantum solitons to resemble closely their classical 
ascendants studied below. 

The equations of motion on the other hand read 

ii> = --~D 2 ^ + A ^ + -k 2 ^*^ - 1)* 
2 4 



jEi = tijdjB - Ji (2.7) 



The Gauss constraint and the remaining Maxwell identity 



- diEi = p B = -eijdiEj (2. 



complete the set of field equations of the model. The charge and current densities 
p and Ji respectively, are given by 

p = $*^_l J i = —[^*D^ - (D^)*^] (2.9) 

2i 

The corresponding energy functional is the sum of four positive terms W = Wd + 
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W b + W e + W v with 



W d = \ J d 2 x\D^\ 2 W b = l - J d 2 x B 2 



^ J d 2 x E 2 W v = ^k 2 J d 2 x (W - l) 2 (2.10) 



W e 

2A 



We will discuss exclusively smooth, particle-like, localized configurations with 
finite energy. Notice that even in the absence of the potential V finiteness of W e 
already implies that these configurations must be neutral. A localized configuration 
with non-zero net electric charge leads to an electric field behaving like 1/r at 
large distances and this makes W e diverge. This neutrality requirement translates 
into a boundary condition for at infinity. In fact |^| must tend to one, since 
otherwise it will lead to infinite charge at large distances. It is possible of course 
to interpolate to a neutral singular configuration but only at an infinite energy 
cost. One could in principle contemplate a logarithmically divergent energy due 
to a non-zero finite charge but it is impossible to allow for an asymptotic value 
of different from one. We conclude that the electrostatic term exactly like the 
standard phenomenological potential term, when combined with the requirement of 
finite energy, leads to the same boundary condition for at infinity. Furthermore, 
we will explicitly demonstrate in the next section that even without V the field 
equations do support the existence of non-trivial soliton solutions. We thus have 
to impose the conditions 

Q = / d 2 x p = (2.11) 



and 

|tf(x)|-»l and |Di*|->0 as |x| -> oo (2.12) 

Such configurations are known [nl to be classified according to the first homotopy 
group of S* 1 into disjoint topological sectors (equivalence classes) characterized by 
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an integer topological charge N, computed from 



N = J~j d 2 x e kl d k ^d^ (2.13) 

For field configurations satisfying the above boundary conditions one may equiva- 
lently use the manifestly gauge-invariant formulas 

N = — [ d 2 xB 

or 

N = ^-j d2x hl(D k ^T(D^) - iB(**V - 1)] (2.14) 
as the integral of the explicitly gauge-invariant topological density 

r = ^fobiOWOW - - 1)] (2.15) 

An argument which distinguishes these three definitions and favours the use of the 
last version will be given in the next section. 

A final technical remark is related to the gauge invariance mentioned above. 
The action (2.4) is invariant under the space dependent gauge transformations 

= exp(iA(-x))^ , A/ = Ai + d{A, Aq = A for arbitrary function A(x). For 
small gauge transformations i.e. for A such that eijdidjA = the topological charge 

is also invariant. Large, and for the U(l) case at hand necessarily singular, 
transformations for which e^didjA = 27rn5(x), A^ changes by n units. To solve the 
field equations one has to eliminate the redundant degrees of freedom by imposing 
a gauge-fixing condition. We choose to work in the gauge V • A = 0. 

Equations (2.7) and (2.8) admit the one parameter family of trivial vacuum 
solutions 

* = exp{ia) Ai = A = (2.16) 

for arbitrary constant value of a. Although the small oscillation analysis around 
a trivial vacuum is of interest by itself and deserves detailed study, the purpose 
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of this work is to analyse the physics of flux-vortices to which we now turn our 
attention. 



3. Static properties 

3.1. Isolated vortices 

We will search for static axially symmetric solutions of the field equations in 
any given non-trivial topological sector. The most general ansatz for such a vortex 
configuration with topological charge N is in polar coordinates (r, 0) given by 

*(x) = X(r)e iN ^ A Q (x) = f(r) A(x) = ^(1 - a(r))4> (3.1) 

At infinity (2.12) implies the boundary conditions 

X{r) -> 1 a{r) -> (3.2) 

Within the above ansatz the electric and the magnetic fields are 

No' 

E = -f'f B = (3.3) 

while the electric charge and current densities take the form 

p = X 2 - 1 J = — «X 2 (3.4) 

r 

Finally the field equations reduce to the following set of ordinary nonlinear differ- 
ential equations 

1 A/" 2 1 

X" + -X' - \xa 2 = -k 2 (X 2 - l)X + 2fX 
r r l 2 

\ if" + If) = -P »" - \a> = aX 2 (3.5) 

The behaviour of the solutions around the origin r = is dictated by the 
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requirement of smoothness and by the equations themselves to be 

X(r)^cir |iV| a(r)^l+c 2 r 2 f \r) -> c 3 + r 2 (3.6) 

where c\, c 2 and C3 are undetermined constants. 

Being unable to solve analytically the equations at hand we proceeded nu- 
merically.After some experimentation we decided that among the various methods 
available the variational approach is the most convenient. The method is based 
on the fact that the field equations (2.7) and (2.8) for static solutions are iden- 
tical to the conditions for the minimization of the energy functional (2.10) under 
the Gauss-law and charge neutrality constraints. This is true even after we insert 
into the energy functional the spherically symmetric ansatz (3.1) [nl and leads in 
that case directly to equations (3.5). We can thus use the variational method and 
minimize the energy in order to find approximate solutions of our equations. 

We approximate the fields of the ansatz by a set of trial functions depending 
on some number of variational parameters. We then evaluate the energy as a 
function of these parameters and minimize the resulting expression. The position 
of the minimum in the parameter space determines the approximate solution, whose 
energy is the value of the energy functional at the minimum. The variational ansatz 
we used, compatible with the boundary conditions at large r is : 

m 1 mi 

X(r) = 1 + e~ kr a(r) = e~ r (3.7) 

;=o ' 1=0 

m 1 , 2m 1 

f'(r) = e-** [£± + Y.c, r -] + e-» fci + J>1] (3.8) 

1=0 ' 1=0 

The asymptotic analysis of (3.5) requires that k — k in the above equations. 
Numerically though it is more accurate to leave k as an independent variational 
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parameter. The optimum value of k is not k. This is no contradiction. The 
asymptotic study actually refers to very large distances, where all fields have al- 
ready reached their vacuum values and do not contribute to the energy. One may 
in principle put k = k but this frustrates the other variational parameters, which 
are forced to take unaturally diverse values in order to lead to the right solution. 
This induces in general higher numerical errors. 

Plugging the above expressions for X and /' into Gauss' constraint equation 
(2.8a) we determine the coefficients c_i, d-\, C{ and d{ in terms of the Z{. The 
boundary conditions at the origin give cto = ol\ — 1, zq = —1 for the N — 1 
vortex. For higher- N vortices one has more conditions due to the faster fall-off of 
X at r = 0. Thus, for the N = 2 vortex we obtain in addition to the above that 
z\ = —k. Finally, the vanishing of the total electric charge eliminates one more of 
the unknown coefficients of the variational ansatz. In our computations as it will 
be explained below, we achieve satisfactory accuracy for m = 10, i.e. with 16 and 
15 independent variational parameters for N — 1 and N = 2 respectively. 

Having solved explicitly Gauss' constraint as well as the neutrality condition, 
we insert the variational ansatz into the energy functional, we compute analytically 
the spatial integrals to end-up with a polynomial of fourth order in the variational 
parameters. A quasi-Newton minimization procedure in then used to evaluate 
its minimum. The method converges fairly rapidly within four or five iterations 
inspite of the large number of parameters. The physical nature of the problem and 
the proper choice of the variational ansatz which restricts the fields to the right 
subspace of the configuration space, are the two factors to which it is reasonable 
to attribute this rapid convergence. 

Figure 1 shows the profiles of the N — 1 vortex solutions for two different 
sets of values of the parameters of the model. In figures la and lb we plot the 
magnitude of the electric current J, the magnetic field B, the charge density p and 
the magnitude of the electric field E devided by A as functions of the radius r from 
the center of the vortex for A = 0.1 and k — 2. Similarly, in figures lc and Id 
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we have the same quantities but for A = 0.1 and k — 0, the case with no quartic 
potential. 

Before we proceed with further results a few comments about the accuracy 
of our numerical computations are in order. To obtain an estimate of the error 
induced by the truncation of the configuration space, we minimized the energy for 
various values of m in order to check the stability of our results against a broadening 
of the trial ansatz (3.7) , (3.8) . The energy at the minimum for m — 10 differed 
from the one for m — 11 in the sixth significant digit. Thus we conclude that 
our energy calculations are correct to within one part in 10 5 or so. To eliminate 
the danger that during a minimization process we were accidentally trapped in a 
local minimum, we repeated each run several times starting with different initial 
configurations. An additional test is provided by the accuracy with which we verify 
the virial relation 

W b = 2W e + W v 

[12] 

obtained by the well known scaling argument of Derrick . We define 

= m. + Wy - m 

and check that in all our calculations A is smaller than 10~ 4 . Finally, the spherical 
symmetry of the ansatz provides another test of the accuracy of our results. Notice 
that the ansatz is spherically symmetric in the sense that a spatial rotation by some 
angle (3 can be compensated by a global internal U(l) rotation by Nf3. This means 
that the corresponding generator Q — NL where L is the field angular momentum 
to be defined in the next section, must vanish for the solution. Q is zero by 
construction and we check, by explicit evaluation of the integral of the angular 
momentum density, that for all solutions L is also zero to within one part in 10 4 . 

Figure 2 shows the values of the ratio R = En = 2/2En=i of the energy of the 
double vortex to twice the energy of the single one, as well as the value of the 
energy of a single vortex itself devided by n for various sets of the parameters 
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of the model. For A = and for static configurations our model reduces to the 
Ginzburg-Landau equations for a superconductor and k corresponds exactly to the 
ratio of the penetration depth to the coherence lenght defined there. 

Figure 2d is in perfect agreement with the results obtained previously in the 

[131 [71 [81 

context of the Ginzburg-Landau theory ' . We divided all energies by n to 
make easier contact with previous results concerning the Ginzburg-Landau theory, 
in which for k — 1 the N = n vortex solution has energy E n = n. The fact 
that the energy of a vortex is an increasing function of the parameters k and A 
is easily understood by the following argument: (we give the argument for k). 
Start with the solution of the static equations for some value of the parameter k 
and a given A which we keep fixed. By construction it satisfies Gauss' constraint 
and is neutral. It can thus be used as a trial initial configuration for the energy 
minimization for a different value k' < k. The energy E(k') of the solution for 
k' is smaller than the value of the energy functional for the above configuration. 
The latter is equal to the energy of the solution for the value k minus the positive 
quantity (1/8) (k 2 - k' 2 ) j d 2 x (^* - l) 2 , where * is the solution for k. Thus, 
E(k') < E(k). Q.E.D. A similar argument with a little extra care due to Gauss' 
constraint applies to A. The other extreme, the no-potential = is shown 

on figure 2c. The critical value of A for the N = 2 vortex to have the same energy 
with two well-separated single vortices is 0.0167. 

The critical line separating the parameter space into regions-I and -II, according 
to the value of R being smaller or larger than one respectively, is shown on figure 
3. Regions-I and -II correspond in the Ginzburg-Landau case to ordinary type-I 
and type-II superconductors. It is interesting to compare the energy of a triple 
vortex to that of three single vortices far from each other. In contrast to the 
Ginzburg-Landau case there exists inside region-I another critical line separating 
the parameter space into two subregions according to whether E(N = 3) is larger 
or smaller than 3 times E(N = 1). Consistent with the well known facts about 
the Ginzburg-Landau model this line too passes through the point (k — 1, A = 0). 
The subregion with E(N = 3) < 3E(N = 1) is further split into two subsubregions 
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depending on the value of the ratio E(N = 4)/4E(N = 1) being larger or smaller 
than one, and so on. 

In the context of a relativistic model we would talk about stable and unstable 
higher- N vortices. For instance, a vortex with N = 2 and energy larger than twice 
the energy of a single one could not be absolutely stable against decay to two 
single vortices. In the model at hand though as we will see in the next chapter a 
completely different picture arises' 11 . 

3.2. The two- vortex system 

Our main goal in this section is to compute the energy of two single vortices 
as a function of their separation d. Again we proceed numerically and choose a 
variational ansatz with the following characteristics: 1) The complex field vanishes 
at the points (±|, 0) on the x-axis. 2) For large d the configuration reduces to two 
well-separated axially symmetric single vortices at a distance d from each other, 
while 3) for d — > it takes the form of an axially symmetric N = 2 vortex located 
at the origin. For the phase O and the magnitude X of the complex field \l/ we 
write [8) 

Q(x,y) = ton" 1 ! 2 2 ? A (3-10) 

x y 4 

X(x) = u X«(|x-^|) X«(|x + i)+ (1-c) (|X ~^ I) 2 (|X+ ^ I) *< 2 >(r) + 5X 

(3.11) 

in terms of the previously determined solutions X* 1 ) and X^> for the N — 1 and 
N = 2 vortices respectively. The function 5X is written in the form 

j j l m 

8X = (|x- -|) (|x+ -|) co S h-\kr) E z ™ ^ ««(2n0) (3.12) 

m=0 n—Q 

uj is a variational parameter, the relative weight of and X^> in the config- 
uration. To fix its optimal value we minimized the energy with respect to uj for 
5X = 0. 
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The variational ansatz is completed by a trial set of functions for the gauge 
potential 

A(x) = coA«(|x-^|) + cuA«(|x + ^|) + (1-cu) A^(r) + <SA(x) (3.13) 
where 5 A = 5A r f + SA^ (f> with 5A r and SA^ defined by 

l m 

5A r = r cosh~ l {r) ^ Yl a ™ ^ sin ( 2n< P) ( 3 - 14 ) 

m=0 n—Q 
I m 

5A^ = rcosh-\r) a ™ r<2m cos{2n<j)) (3.15) 

m—0 n—Q 

As mentioned above d = (d, 0) to place the two vortices at (±|, 0) along the x 
axis. The angular dependence of the fields was restricted to the form given above 
by the requirement on the configuration to be invariant under spatial reflection 
with respect to either one of the coordinate axes. 

We used the same numerical procedure as in the isolated vortex case to de- 
termine the optimum values of the remaining variational parameters afj,a^j, Zij 
and k. The dominant source of error in our results is due to the truncation of the 
configuration space. Again an estimate is obtained from the change in the value 
of the quantity of interest as we vary the number of variational parameters. To 
achieve an error of the order of 0.1% in the total energy for all values of d, it was 
enough to set I — 1 for d > 3, while for d < 3 it was necessary to take I = 2 i.e. 
a richer variational ansatz with 18 parameters. The uncertainty in W e was a little 
larger, something like 0.2 — 0.3%, but since W e contributes always a small fraction 
to the total energy, this error is negligible. Contrary to the isolated vortex case 
we did not perform the spatial integrations analytically. Instead we carried out 
these integrations numerically. Appropriate choice of the grid and the boundaries 
reduced the corresponding error to the order of 0.01%. 
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In figure 4 we plot the energy divided by tx of the two-vortex system as a 
function of their separation d for four sets of the parameters k and A correspond- 
ing to the points A, B, C and D shown in figure 3, from well inside region-II, to 
somewhere in region-I. We have included the point C on the critical line with 
(k — 0.5 , A = 0.0115). Figure 4a corresponds to parameter values far from the 
critical line in region-II. The force between the vortices is everywhere repulsive. 
In all other cases one immediately recognizes a region of attraction and a region 
of repulsion of the two solitons. The local extremum in the energy leads to the 
exciting possibility of the existence of dumb-bell shaped doubly-charged static so- 
lutions of the field equations for a rather wide range of parameters. We will have 
the opportunity to elaborate on this in a future publication 1141 . 

Finally, figures 5, 6, 7 and 8 show the total energy density wt, the gauge- 
invariant topological density r, the magnetic field B and the electrostatic energy 
distribution w e respectively, of the optimum two-vortex configuration correspond- 
ing to the parameter values (« = 0.5 , A = 0.0115) on the critical line. The 
pictures in each one of them refer to distances d — 0, 2, 4, 6. Notice that of the 
two definitions B/2n and r of the topological density the former follows less closely 
the energy distribution and as such it offers a bad description of the two vortex 
positions. We have checked, although we do not show it here, that the topological 
density in terms of the complex field alone is even worse in that respect. Notice also 
that the distribution of w e does not follow that of the total energy. It spreads over 
the whole region between the two vortices and is actually responsible for the bump 
in their interaction energy. Being so small in magnitude it does not alter the final 
picture that the two vortices can be safely considered far from each other already 
for d = 4. These features of w e as well as its negligible magnitude compared to w t 
are also shown on figure 9 for the case of an axially symmetric N — 1 configuration 
with k = 0.5 and A = 0.0115. The virial relation given in the previous section 
shows that W e is never dominant. It can at most be 1/3 of Wt — for any value 
of the parameters. 

We end this section with a remark about the numerical procedure used to 
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compute W e at each step of the iterative minimization process. We exploited the 
linearity of the Poisson equation to compute once and for all the electric field 
created by the charge density due to the configuration (3.10) to (3.12) in terms of 
the parameters of the ansatz. Thus the calculation of the electric field after each 
iteration is immediate. One does not have to solve the Poisson equation again but 
just plug into the general formula the new values of the parameters. W e is finally 
given by a straightforward spatial integration. This way the computation is more 
accurate and much faster. 



4. Vortex Dynamics 
4.1. Canonical structure of the model 

In order to investigate the dynamics of vortices we start with the canonical 

[15] 

structure of the model. We follow the standard procedure ' applied to gauge the- 
ories to determine its fundamental Poisson brackets. As mentioned in a previous 
section to eliminate the redundant degrees of freedom due to the gauge invariance 
of the model we impose the condition 

V- A = 



Since the action does not depend on Aq, Aq is a dependent variable satisfying 
the Gauss law constraint which we solve to obtain 

A (x, t) = -— [ d 2 x' ln\x - x'| (^(x ; , f) - 1) (4.1) 
2tt J 

(use was made of V 2 ln\x. — x'| = 27r5(x — x')). With the above gauge-fixing 
condition the purely electromagnetic part of the action takes the form 

S A = l -Jd 2 x[ i(E r 2 + E L 2 ) - B 2 ] (4.2) 

where the transverse and the longitudinal parts of the electric field are E r = —A 
and Ex, = — VAq respectively. 
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The canonical momentum Wk conjugate to is then ir^ = dC/dA^ = —ET/X, 
while \l/ and n = dC/d^ = i^fr* are conjugate of each other. Thus, the fundamental 
Poisson brackets are 

{^(x,t),M/*(x / ,t)} = -^(x-x / ) 
{Aii^t),^'^)} = A5 T y(x-x') (4.3) 

{tt(x,t), A(x',i)} = = {*(x,t),i i (x',t)} 

where the transverse 5 function is defined as 5 T y(x— x') = — V _2 <9i<9j) 5(x— x'). 
Finally, the Hamiltonian (energy) of the system is 

H = J d 2 x[ i-(E T 2 + E L 2 ) + l -B 2 + i|D^| 2 + i/t 2 (^ - l) 2 ] (4.4) 

All other Poisson brackets are obtained from the ones given above. It is 
straightforward to check that 

{A (x,t),^(x / ,t)} = -^/nlx-x'l tf(x',t) 

{A)(x,t), A t (x,t)} = = {A (x,t), ii(x',t)} (4.5) 

as well as that the equations of motion (2.7) and (2.8) coincide with Hamilton's 
equations 

*(x,t) = {tf(x,t),tf} 
A fc (x,0 = {A fc (x,*), J ff} (4.6) 

7Tjfe(x,t) = {lT k (x,t),H} 

We are now ready to proceed with the construction of the field momentum 
and angular momentum for the model at hand, valid in all topological sectors. We 
will show' 11 191 that the corresponding conservation laws describe the basic feature of 
vortex motion in the absence of external forces, namely spontaneous pinning. 
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4.2. Momentum and angular momentum in all sectors 



The Noether expression for the linear momentum in our model is 




d 2 x (-7rd k ty - TTjd k Aj) 



(4.7) 



It is derived under the assumption that all fields approach their vacuum values fast 
enough as r — > oo so that all integrals are meaningfull and surface terms appearing 
in intermediate steps are zero. It is straightforward to check that the above ex- 
pression is ill-defined for any vortex configuration with non-zero topological charge. 
Under the same assumptions though which led to (4.7) the latter can be brought 
to the form 



which is well-defined for any smooth finite-energy field configuration with arbitrary 
topological charge. The two expressions (4.7) and (4.8) differ by exactly those 
surface terms which were omitted from the former and make the latter finite and 
unambiguous in all topological sectors. 

All the defining properties of the momentum are verified for P k . This is no 
surprise. They are formally valid for Pj? and this differs from P k only by surface 
terms. In any case, one can directly and unambiguously show using the equations 
of motion (2.7), (2.8) and the Poisson brackets derived above first, that P k is 
conserved 

s^ = ° (49) 

second, that it is the generator of spatial displacements 



and finally, that it is gauge-invariant. The only condition on the configuration 
T = (^, Aq, Ai) necessary for the derivation is that Aq vanishes faster than r _1 





{P k ,T} = d k T 



(4.10) 
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as r — > oo. All our vortices have Aq approaching zero at infinity exponentially fast 
and consequently they safely belong to this set of configurations. Furthermore, 
using Gauss' constraint we can rewrite the momentum in the manifestly gauge- 
invariant form 111 

Pk = e ki J d 2 x (2tt Xi r + i EiB) (4.11) 

The second term is the Poynting vector, the pure electromagnetic contribution to 
the field momentum. 

Thus, independently of derivation, P k is the correct form of the momentum 
in our theory, valid in any topological sector and reducing to the naive expression 
in the trivial sector N = and for configurations allowing for free integrations by 
parts. 

Before we discuss the implications of the above formula of the linear momentum 
on the motion of vortices, we would like to construct in a similar way the correct 
form of the angular momentum in our model. Again, the Noether expression 

L N = - J d 2 x [ Xi (irdj^ + n k djA k ) + mAj ] (4.12) 

is formally conserved and generates rotations but is divergent in any non-trivial 
sector. Allow for integrations by parts and ignore surface terms to rewrite it in the 
form 

L — J d 2 x[ ix 2 {di-Kdfi + diTx k djA k ) + ] (4.13) 

or as a manifestly gauge invariant quantity in terms of the topological density r 
defined in section 2 

L = - J d 2 x (ttxV + i x • EB) (4.14) 

Using either one of the last two formulas one can show that for essentially 
all finite-energy configurations as in the case of the linear momentum, L is well- 
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defined, conserved and generates spatial rotations, i.e. it satisfies 



{L,tf(x,t)} = djXidjV 

{L, Ajfe(x, t)} = eijXidjAk + e ki A x (4.15) 

L is then the proper definition of the angular momentum when dealing with 
topologically non-trivial configurations. As already explained it vanishes for the 
spherically symmetric neutral solutions studied in section three. 

4.3. Vortex motion 

In a following paper 1141 we will present details of a numerical simulation of the 
motion of the flux vortices described above under the influence of various external 
forces. Here we would like for the sake of completeness to predict analytically and 
without any approximation the essential features of their dynamical behaviour' 11 . It 
is already apparent that the momentum defined in the previous section is not a 
measure of the translational motion of a vortex but instead it describes its position. 
In fact the momentum of a static axially symmetric vortex with charge N centered 
at a is Pk = 27iNeki(ii. It is also clear that a localized free vortex of arbitrary shape 
moving in formation at constant velocity V{ would have = 27r7Vefcj(a° + Vit). For 
any v i ^ this is forbidden by the linear momentum conservation law. A free 
vortex is spontaneously pinned as a consequence of momentum conservation. 

Let us define the "guiding center" R of a generic configuration with N ^ by 

Ri = ~^ijPj (4-16) 

Since under a rigid displacement of the whole configuration by c, R changes to 
R + c and for a spherical vortex it coincides with its geometric center, we naturally 
interpret it as the mean position of the configuration. Note also in support of this 
interpretation, that R is related to the first moment of the topological density r, 
so that for nearly spherical configurations R is close to their "center of topology". 
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a. In the absence of external forces both N and P are conserved and 

^-R=0 (4.17) 
dt v ; 

A generic vortex-like configuration produced in the system will of course fluctuate 
in its details but it will remain pinned at its initial mean position. 

b. Consider next the response of such a vortex to an externally prescribed 
electric current Ji(x, t). Its effect is studied by adding to the action the term 
SS = J d 2 x dt Ji(x, t)Ai. For concistency of the model the external current should 
have zero divergence. The new term in the action modifies the A{ equation of 
motion (2.7) by the substitution Jj — > Jj + Jj on the right-hand side. Because 
of the external current the momentum and the angular momentum are no longer 
conserved. Instead, their time derivatives are equal to the external force and torque 
respectively. Indeed one can easily verify Newton's equation: 



d_ 
~dt 



P k = Ft orentz = - J d 2 x e kl J,(x, t)B(x, t) (4.18) 



as well as 



-L = Torque = J d 2 x xJiB (4.19) 



d 

~di 

Naturally the force on the vortex is opposite to the Lorentz force acting on the 
external current, while the torque is by definition expressed in terms of the force 
density f k = —e k iJ[B as the integral of tijXifj . We now use (4.16) and the fact 
that for arbitrary J N is still conserved to obtain 

dR k 1 



dt ' 2nN 



d z x J k B (4.20) 



In the idealized situation of a homogeneous throughout the vortex external 
current the above formula simplifies to 

dRk 



dt 



-■hit) (4.21) 



In general, the naive expectation based on the usual Newtonian reasoning and 
Galilean invariance would be that the vortex should accelerate in the direction of 
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the force acting on it. Instead, (4.16) combined with dP^/dt = for a generic 
force F}~, shows that the equation of motion of the mean position of the vortex is 

^-R k = -— !— e kl Fi (4.22) 
dt 2nN y ' 

i.e. the vortex moves with speed |F|/27r|iV| and at ± 90° relative to the force for 
positive or negative N respectively. In the special case of the homogeneous external 
current considered here, the force is also proportional to iV and all vortices move 
with the same velocity equal to minus the external current itself. If in particular the 
latter is due entirely to the condensate charges taken here by convention positive 
with unit charge density, its value is exactly equal to minus the velocity of the 
carriers. Thus the vortex will reorganize itself during a transient period following 
the onset of the external current and it will move with the same speed but opposite 
to the current carriers. Furthermore it should be pointed-out that in the context of 
our field theory model no approximation other than the implicit assumption that 
the vortex remains localized was necessary. The position interpretation (4.16) of 
the momentum converts Newton's law into an equation giving directly the velocity 
of the vortex in terms of the applied force, while the eij of (4.16) makes the vortex 
move at an angle (N/\N\) 90° relative to the direction of the external force (Hall- 
behaviour) 111 [9] . The above general conclusion, reached without ever solving an 
equation of motion, does not depend on the details of the Hamiltonian of the 
system. Any potential V or any additional higher spatial derivative terms in H 
modifies the detailed profile of the vortex solutions but it does not alter their 
dynamics. 

The guiding center interpretation of the momentum (4.16) is not applicable in 
the case of topologically trivial (N = 0) configurations. Newton's law (4.18) is of 
course still true for any external force, but this does not strictly speaking tell us 
much about the actual motion of the vortex. In contrast to the case of ordinary 
particle dynamics with Poincare or even Galilean invariance the relation between 
momentum and velocity is not a priori known in the present model. 
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5. Discussion 



We presented a general treatment of the motion of flux-vortices in a two pa- 
rameter field theory model describing as stated in the introduction a rather wide 
class of idealized physical systems. We argued that independently of the details 
of the Hamiltonian the vortices exhibit Hall-behaviour, like electrons moving on a 
plane under the simultaneous action of a perpendicular magnetic and an in-plane 
electric field. The picture is identical to the one derived previously for the dy- 
namics of magnetic bubbles in the ferromagnetic continuum' 91 . The fact that the 
solitons in these two field theories are topologically different, classified by the sec- 
ond homotopy group of S* 2 in one case and by the first homotopy group of S* 1 in 
the other, does not seem to play any role in our discussion. 

An external current was shown from first principles to pull the vortex in the 
opposite direction. It is not clear at this point what is the relation of our model 
with the ones considered previously' 41 [16) in connection with superconductivity. It 
is tempting to think of it as a field- theoretic realization of the hydrodynamic 
model discussed in reference [16] in the vanishing friction limit. The conclusion 
though reached there about vortex motion, based on the phenomenological use of 

[17] [18] ~* 

the Magnus force known from fluid mechanics , is V = J. This is opposite to 
ours, which is correspondingly reminiscent of the opposite sign Hall effect reported 

[3] 

in high-T c superconductors 

The Hall motion of the vortices derived above is attributed to the radical 
change in the role of the momentum of the theory as a result of the underlying 
topology. The Hall behaviour of an isolated vortex is exactly due to the fact that 
the linear momentum (4.11) contains a piece which is equal to the first moment of 
the topological density. An immediate consequence of this fact combined with the 
fundamental property (4.10) of the momentum is that 

{P h P 2 } = 2irN (5.1) 
We thus conclude that necessary condition for the manifestation of Hall motion 
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is that the translation part of the symmetry algebra of the model admits central 
extension. It is well-known [191 that this is not possible in the case of the Galilean 
or the Poincare algebra in any spatial dimension higher than or equal to two. It is 
not even true for the Euclidean algebra E(D) for D > 3. Thus, the Hall behaviour 
encountered above is expected a priori only in two- dimensional systems with spa- 
tial £'(2) algebra, or in D-dimensional systems with translational symmetry alone. 
In the system at hand one can immediately check that 

{L,P 1 } = P 2 {L,P 2 } = -P X (5.2) 

and verify that it belongs to the first category as expected. 

We would also like to stress that although (4.8) is true even in theories with 
canonical momentum 7r proportional to the time derivative of the complex field, our 
previous reasoning does not go through since P is then related to the first moment 
not of the topological density but of a quantity which actually vanishes for static 
configurations. Thus, another condition which is necessary is that the equations of 
motion are first order in the fields which carry the topology in the particular model. 
In other words, we do not expect Hall- motion of the vortices in the relativistic-like 
model 1201 derived as a dynamical extension of the Ginzburg-Landau theory of a 
superconductor appropriate for temperatures far below T c . 

On the other hand, the above formalism applies to the motion of vortices 
in a generalization of the present model which allows for dynamical positive-ions 
forming a lattice, of the (N^, N ss ) with N ss ^ — 2Nf l vortices of the high-T c super- 
conductor model proposed in reference [21], as well as to the motion of vortices in 

[22] 

easy-plane ferromagnetic films , to mention just a few systems of considerable 
interest. 
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FIGURE CAPTIONS 

Figure 1: N — 1 vortex profiles, (la) and (lb) correspond to k = 2 and 
A = 0.1, while (lc) and (Id) to k = and A = 0.1 

Figure 2: Plots of R = En=2/2En = i and of En=i/it as functions of k or A for 
the values of the second parameter as shown. (2d) corresponds to the Ginzburg- 
Landau model, while (2c) to the no- potential case. 

Figure 3: Regions-I and -II of the parameter space. 

Figure 4: The energy of the two- vortex system as a function of their separation 
for the values of parameters corresponding to the points A, B, C and D of figure 
3. The dashed line is drawn at twice the energy of a single vortex. 

Figure 5: The total energy density of the two vortices d = 0, 2, 4 and 6. 
k = 0.5 and A = 0.0115. 

Figure 6: The topological charge density r of the two vortices for d — 0, 2, 4 
and 6. ac — 0.5 and A = 0.0115. 

Figure 7: The magnetic field of the two vortices for d — 0, 2, 4 and 6. ac = 0.5 
and A = 0.0115. 

Figure 8: The electrostatic energy density of the two vortices d — 0, 2, 4 and 
6. ac = 0.5 and A = 0.0115. 

Figure 9: Plot of the densities w t and 150u> e for an isolated N — 1 vortex. 
k = 0.5 and A = 0.0115. 
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